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Parker's analytic Cartesian interface dynamo is generalized to the case of a shear layer of finite thickness and low resistivity ( "tachocline" ) , 
bounded by a perfect conductor ("radiative zone") on the one side, and by a highly diffusive medium ("convective zone") supporting 
an a-effect on the other side. In the limit of high diffusivity contrast between the shear layer and the diffusive medium, thought to be 
relevant for the Sun, a pair of exact dispersion relations for the growth rate and frequency of dynamo modes is analytically derived. 
Graphic solution of the dispersion relations displays a somewhat unexpected, non-monotonic behaviour, the mathematical origin of which 
is elucidated. The dependence of the results on the parameter values (dynamo number and shear layer thickness) is investigated. The 
implications of this result for the solar dynamo problem are discussed. 
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1 Introduction 

Interface dynamos are a widely discussed class of astrophysical dynamos, especially in the solar context 



(Petrovay 2000 Solanki et al. 2006, Charbonneau 2005). In contrast to the textbook case of Parker's 



migratory dynamo, in an interface dynamo the a and Q effects operate in spatially distinct but adjacent 
layers. In this setup the dynamo wave will arise as a surface wave propagating along the interface of the 
two regions. While dynamo problems with spatially separated a and O had been considered since the 



1970's (see e.g. the references in Zeldovich et al. 1984), the classic analytic formulation of the problem is 



dueto|Parkerj([1993|). 

It has long been suspected ( Parker || 19 75 



numerical simulations ( Browning et al. 2006 



Schiissler 1984, Petrovay 1991), and recently demonstrated in 



that the toroidal magnetic flux generated in the solar dynamo 
can only be stored below the base of the convective zone proper, as various flux transport effects (buoyancy 
and pumping) remove it quite effectively from the convective zone. With the discovery of the tachocline 
layer it became clear that this hypothetical flux reservoir coincides with the strongest rotational shear in 
the solar interior, which is plausibly also responsible for the production of the strong toroidal magnetic 
fields, i.e. that this is the site of the O-effect. The site and nature of the a-effect is much less clear; in 
interface dynamo models it is assumed to be concentrated in the deep convective zone, just above the 
tachocline. 

Current helioseismic evidence indicates that at low heliographic latitudes the tachocline is situated 
immediately below the adiabatically stratified convective zone. (At higher latitudes some overlap may be 
present.) This suggests that the tachocline is much less turbulent than the convective zone proper. This is 
also consistent with another physical consideration: the high overall toroidal field strength that must be 
present if all the magnetic flux emerging during a solar cycle resides in the thin tachocline layer should 
strongly suppress turbulence there. The strong subadiabatic stratification of the tachocline will then inhibit 
the penetration of meridional circulation to it, presenting a difficulty for flux transport models, the main 
contender of interface dynamo models. 
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There is some disagreement over whether the oscillatory magnetic field generated by the dynamo pervades 
the whole of the tachocline or not. For effective magnetic diffusivities less than 10 8 cm 2 /s or so, the pene- 
tration depth (skin depth) of the oscillatory magnetic field is much less than the helioseismically inferred 
tachocline thickness (order of 10 Mm), so the penetration cannot be complete. In this "slow tachocline" 
scenario ( Garaud|2001 Brun and Zahn|2006 ) the dynamical changes in the bulk of the tachocline take place 
on the long diffusive timescale. In the alternative "fast tachocline" scenario (Forgacs-Dajka and Petrovay 



2001, 2002, Forgacs-Dajka 2003), hydrodynamical instabilities maintain some modest level of turbulence 



in the tachocline, so the resulting higher skin depth allows the dynamo field to penetrate deeper. In this 
case it is the Lorentz force in the dynamo generated field that limits the penetration of differential rota- 
tion deeper into the radiative interior, so tachocline thickness and skin depth are intimately related and 



should agree within an order of unity factor ( Petrovay 2003 ) . In the present paper we will consider the 



fast tachocline case; some observational support for this scenario comes from the recent detection of solar 
cycle related changes in tachocline properties ( Baldner and Basu]|2008 ) , as predicted by the fast tachocline 
models. 

The interface dynamo option for the solar dynamo is attractive, as it would explain the high (~ 100 kG) 
field strength of toroidal fields and their equatorward propagation, while still being consistent with phys- 
ical assumptions about the magnetic diffusivity values in both the convective zone and tachocline and 
about the meridional flow amplitude. Interface dynamo models for the Sun have been constructed, e.g., 



by Charbonneau and MacGregor (1997), Tobias (1996) or Markiel and Thomas (1999). Unfortunately 



none of these models show a compelling detailed agreement with observed features of the cycle. Their 
main alternative, the flux transport dynamos ( Dikpati and Charbonneau|1999 Chatterjee et a/.j|2004 ), in 
contrast, are more straightforward to parameterize so that they can reproduce some salient features of the 
solar cycle reasonably well; however, the physical consistency of the models is very much in doubt, espe- 



cially regarding their assumptions on diffusive and advective magnetic flux transport processes ( Choudhuri 



2008). 



The now classic analytic Parker (1993) interface dynamo consisted of two adjacent semiinfinite domains: 



one highly diffusive domain home to the a-effect and another, low diffusivity domain characterized by 
strong shear. As we have seen, the solar tachocline, corresponding to the shear layer in the Cartesian 
interface dynamo, has a quite limited thickness, so it is not particularly well represented by a semiinfinite 
domain. The aim of the present paper is to generalize Parker's analytical interface dynamo to the case 
when the shear layer has a finite thickness. 



What makes this problem particularly relevant is the recent proposal (Petrovay 2007) that combining 



an interface model with a fast tachocline scenario both the magnetic field amplitude and the thickness of 
the tachocline can be determined (in the simplest case, by solving two algebraic equations). Under some 
plausible assumptions, two solutions are found: one with strong field and thin tachocline and another 
one with weaker field and a thicker tachocline. These represent "high" and "low" states of solar activity, 
suggesting a link with the phenomenon of grand minima. Based on an analogy with surface gravity waves, 
this model assumed that a finite shear layer thickness h reduces the growth rate of the interface dynamo 
by a factor tanh(hk), where k is the wave number. One objective of the present work is to examine the 
validity of this assumption. 

We note that interface dynamos with finite shear layers have previously been considered in many nu- 
merical implementations (e.g. Charbonneau and MacGregor 1997, Markiel and Thomas 1999) as well as 



in a semianalytic Cartesian setup more general than our current model (Zhang et al. 2004). The effect 
of varying h, however, was not studied in any of those models, due to the large number of other free 
parameters and to the different focus of those papers. It is therefore legitimate to consider a setup where 
h is a variable parameter while in other respects the model is kept as simple as possible. 

The structure of our paper is as follows. In section 2 we analytically derive the dispersion relations in the 
limit relevant to the Sun. Section 3 presents graphical solutions of these dispersion relations and discusses 
their properties. Finally section 4 concludes the paper. 
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2 Derivation of the Dispersion Relations 

2.1 Generic complex dispersion relations 

We use a Cartesian setup where the coordinate z corresponds to height above the bottom of the convective 
zone in the solar application, while x and y correspond to heliographic latitude and longitude, respectively. 



Let r\ denote the magnetic diffusivity, a the a parameter of dynamo theory (see e.g. Petrovay 2000) and 
$7 = dvy/dz the shear rate of the velocity field representing radial differential rotation. The magnetic field 
is decomposed into a "toroidal" component B y and a "poloidal" component V x A y e y where A y is the 
toroidal vector potential. 

Consider now the following 3-layer "sandwich" setup: 

Convective zone z > r\ = r] + q^0 O = B y = B(x, z, t) A y = A(x, z, t) 

Tachocline — h < z < r\ = r\_ q = f2^0 B y = b(x, z, t) A y = a(x, z, t) 

Radiative interior z < —h ry = a = f2 = B y = A y = 

The parameters r/ + , r/_, a and f2 are assumed to be constant in their respective layers. 
The a£l dynamo equations in z > then read 

OB fd 2 B 8 2 B\ n ,., , 

+ ^H = 0, (la) 



dt \ dx 2 dz 2 

dA (d 2 A 8 2 A\ 

+ TOT = ccB, (lb) 



while in the tachocline we have 



dt \ dx 2 dz 2 



db fd 2 b d 2 b\ „da 



dt V - [dx- 2 + dz~ 2 ) = n dx-> <2a) 

da (d 2 a d 2 a 

dt \ dx 2 dz 



V- ( 7SZ2 + 7W2 ) = °- ( 2b ) 



The matching conditions at z = are 



, r, A da dA db dB . , ,. 

6 = 5, = A, ^r = ^-, r ?-^r= r ?+^T> (3a,b,c,d) 



while at z = — h they simply read 



dz dz ' <9z dz 



0, a = 0. (4a,b) 



Following the standard procedure for the analytical solution of interface dynamo equations (Parker 1993 



Petrovay and Kerekes|2004 ) , for the solutions in terms of our variables f(x, z, t) we consider normal modes 



of the form f(z) exp[(<r + iw)t + \kx\. In z > 0, the z-dependent part f(z) is sought in the form 

B = Ce~ Rz A=(D + Ez)e~ Rz , (5a,b) 

where R = S + iQ, and C, a, u, k, S and Q are all real, while D and E are complex. As we expect the 
field to vanish as z — > oo, S must be positive. 

In the shear layer the solutions have a similar form, but, owing to the finite thickness of the layer, the 
solution is in general a superposition of modes growing and decaying with z: 

a = (Je rz + Jie- r "), b = (L + Mz)e rz + (L x + M lZ )e- rz , (6a,b) 
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where r = s + iq with s and q real and s > 0, as before. Equations (|4a,b) then imply 

J 1 = -Je~ 2rh , Li = /iMi + (hM - L)e~ 2hr , (7a,b) 

respectively. 

Substitution of our trial solutions (f5k,b), (J6k,b) into the four interface fitting conditions (pk-d) yields 

J(l - e" 2 ^) = D, (8a) 

L(l - e~ 2hr ) + h(Me~ 2hr + M^ = C, (8b) 

rJ(l + e- 2hr ) = E-DR, (8c) 

rL(l + e~ 2hr ) + M + M 1 - rh(Me~ 2hr + Mi)l = -rj + RC. (8d) 

On the other hand, substituting (7i,b) into the trial solutions (J5^,,b) and ([6^,b), and substituting these 
into the dynamo equations ([jji,b), (2i,b) we obtain the dispersion relations 

r) + R = a + \oj + r/ + k , (9a) 

■q_r =a + iuj + r]-k , (9b) 

E = aC/2r] + R, (9c) 

-27]^r(Me rz - M ie ~ rz ) = ikn,J(e rz - e- rz e~ 2hr ) (9d) 

the last of which implies the two separate relations 

M = -lkQJ/2r]_r, M l = Me~ 2hr . (10a,b) 

Using the relations ([8^-c), as well as ([9j3,d), equation (|8]l) can be written as 

rR(fi 2 r + R5)(r + R5)/k A = \\N[8 - (1 - 5 2 )hr], (11) 



where we introduced the dynamo number N = aVL/rj^k^ and the notation 5 = tanh(/ir). Equation (11) 
establishes a relation between the complex vertical wavenumbers r and R. A second such relation is derived 
by eliminating (a + iw) from equations (|9p,,b): 

R 2 = fi 2 r 2 - Gu 2 - l)fc 2 , (12) 

where the diffusivity contrast /x 2 = r]—/r]+ was introduced. 

Given the parameters N, \i and b (or, equivalently, h) defining the problem, the pair of complex dispersion 
relations ( 11 ), (|12[) can in principle be solved for r and R. The frequency ui and growth rate a of the dynamo 



wave then follow from (|9a,b). The nine parameters C, D, E, J, J\, L, L\, M and M\ are then determined 
by the relations ([7]H([8]:) and (|9p,d). Substituting these into (|5^,b), ([6^,b) together with the derived values 
of r and R finally should yield the full z-dependence of the solution for a given value of k. 

2.2 Solar limit 

Unfortunately, an attempt to actually implement this procedure analytically encounters difficulties. Equa- 



tions (11), (12) look deceivingly simple, yet the actual calculation soon becomes quite involved. To make 
progress, from this point onwards we restrict attention to the limit relevant for the solar interior. 

In the solar convective zone the dynamo wave number and frequency are empirically known to be of 
order k ~ 1/R® and to ~ rj + /R% ) ~ f]+k 2 , where R® is the solar radius. In the tachocline below, the 
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magnetic diffusivity is uncertain but surely several orders magnitude lower than the fixed value 77+: we 
may then safely take the limit r\- — > 0. With the fixed values of r/+, k and ui given above, the penetration 
depth (skin depth) of the oscillatory dynamo magnetic field will be (2r/_/a;) 1 ' 2 ~ fi>/k, which also implies 
\r\ ~ k/fi. As the thickness h of a fast tachocline is expected to be a few skin depths, so in this limit we 
also expect h ~ fi/k — > 0, suggesting the use of a fixed nondimensional layer depth H = hk/fi instead. 



Finally, inspection of equation (11) shows that the leading terms on the left-hand side will be of order 
1//U 2 , while those on the right-hand side are of order N. So in order to keep the dynamo supercritical in 
the solar limit we need to keep fj?N fixed. 

Summing up these considerations, the limit relevant to the solar case is 

/i — > 0, while H = hk/fi = finite and J\f = /i N = finite. (13) 



Equation (12) then simplifies to 

R=( jl 2 r 2 + k 2 ) 1 ' 2 . (14) 



Substituting this into (11) yields an equation for the complex wavenumber r only. In analogy with the 
nondimensional layer depth H, we now introduce the nondimensional equivalent of r as ^ = fir/k: clearly, 
^ will remain finite in the solar limit, as discussed above. Expressing r with ^ and keeping only leading 



order terms in the limit of [i — > 0, equation (11) finally takes the form 



[* 2 (<5 2 - l) 2 H 2 + 2^S(5 2 - 1)H + 5 2 ] Af 2 /16 + i^ 2 (* 2 - 1)5 [$>(5 2 - l)H + 5] Af/2 = 5 2 ^{1 + ^ 2 ) 2 . 

(15) 
In this equation, ^ and 5 are complex, while the parameters Af and H are real. We write 5 and \& explicitly 

as 

5 = 5 R + iSi, ^ 2 = a + iv, (16a,b) 

where 

a = cr/r] + k , v = oj/t] + k (17a,b) 



are the dimensionless growth rate and frequency of the dynamo wave, and (U6p) follows from d9p) in the 
limit [i — > 0. 



Upon performing these substitutions, equation (15) is split into real and imaginary parts. (The calcula- 



tion is very tedious but still manageable by numerical algebra packages.) This finally yields two dispersion 
relations for the nondimensional growth rate a = f3 — 1 and the nondimensional frequency v. These are of 
the form 

f 1 (a,v;H,M) = 0, / 2 (a, u; H,Af) = 0, (18a,b) 

where 

h{a,v\H,N) = (a 11 H 2 -a 12 )M 2 + b 1 M + c 1 , f 2 (a,v;H,N) = (a 21 H 2 - a 22 )M 2 + b 2 N + c 2 (19a,b) 
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an = i/(l + 4 + <5|) - 4*^,(1 - 4 + $) - 2i/(<5£ + 35 2 R 5J - dj), 
ai2 = 2o/6\r, 

6i = 8[2i/(l + 2a)5 R 5! - (a 2 + a - u 2 ){5 2 R - 5])}, 

Cl = 32[(a 2 + a- u 2 )5 R - u(2a + 1)<5 7 ] [{a 2 + a- v 2 )bi + i/(2er + 1)or], 
a 21 = cr(l + 4 + <5f) - 4^(5/(1 - 4 + 5?) - 2a(4 + 3<5|<5| - of), 

c-2 c2 

«22 = o R -dj, 
b 2 = 8[l/(l + 2or)(<5| - 5?) + 2(a 2 + or - ^)Wj], 
C2 = I6[cr(cf + l)(<J fl - (J/) - v{v + l)fo + v{y - \)bi\ 
x [er(or + l)(5 R + <5 7 ) - u(u - l)5 R + 1/(1/ + 1)<J/]. 

Recalling that 5 = tanh(/ir) = taiah(Hfy), the real and imaginary parts of 5 are 

cosh(iJs) sinh(/7s) cos(Hq) sin(Hq) 



Sr 



$i 



cosh 2 (Hs) - sin 2 (Hq) ' ' cosh 2 (iJs) - sin 2 (Hq) 

where s = jis/k and q = jj,q/k are related to a and v by virtue of equations (f9k,b): 



(20a) 
(20b) 
(20c) 
(20d) 
(20e) 
(20f) 
(20g) 

(20h) 



(21a,b) 



,[a + {d 2 + v 2 )y% 



f 



a + [a 2 + u 2 } 1 / 2 ' 



It is straightforward to check that, in the limit H — > oo (i.e. 5 — > 1), Parker's original solution 
(cf. equations (60), (61) in Parker 1993). 



(22a,b) 
is retrieved 



3 Solution 



While the equations (18i,b) are quadratic in J\f, they are in general extremely complicated functions of 
the other variables. Nevertheless, equations (18i,b) lend themselves most readily to a graphic solution. 
Fixing the value of the control parameters N and H, solutions correspond to the crossing points of the 
zero contours of functions f\ and fi in the va plane. Then, varying H while keeping either N or the total 
shear fixed one can study how the solutions deviate from the Parker solution (H — > oo limit) for finite 
values of H . 

For the purposes of this study we take the "textbook" case M = 12V2, already considered by Parker 
(1993). In the limit H — > oo, considered by Parker, there is only one growing mode characterized by a = 1 
and v = v2, as borne out in figure [lp. As H is decreased, the solution bifurcates, and two sequences of 
apparently growing modes appear. Recall, however, that during the derivation of the dispersion relations 
in section 2, several steps may have resulted in the appearance of spurious solutions. Indeed, on physical 
grounds we expect to return Parker's solution in the H — > oo limit but also to behave as a — > at some 
finite positive value of H . On the basis of this consideration we must reject the "upper" (i.e. higher growth 
rate) sequence of solutions in figure [TJ converging to a = 1.67, v = 2.01, as spurious. The lower sequence, 
representing physical growing modes, in turn disappears for H < 1.45, i.e. the dynamo is suppressed for 
such low layer depths. 

The variation of a with H is shown in figure [2} The variation can be roughly represented by a shifted 
tangent hyperbolic function. This was to be expected in analogy with surface gravity waves (Petrovay 2007), 



the shift here being due to the suppression of growing modes at a final layer depth. The detailed behaviour of 
the solution is, however, considerably richer that this simple theoretical expectation, displaying a surprising 
non-monotonic modulation. This is a real physical interference effect, due to the interaction of the waves 



July 12, 2010 0:19 Geophysical and Astrophysical Fluid Dynamics findcpth'rcv2 



Interface dynamo over finite shear layer 





JT= 16.971 


# = 50.00 


2.0 










; / 


1.5 


\ 


' / 


ib 1.0 


/ 
. . . . i . . 


\ """' 


0.5 




^v 


0.0 


^r, , i , , , , i , , , , 


i i 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

v 







JT= 


16.971 




H= 2.00 


2.0 






■ ■ V 1 ■ ■ 


1 r 


/ 


1.5 




^---^ 


-'' 




y 


Zt> 1.0 








+ P 


,--'' 


0.5 








/l 
/ 1 
/ 1 

/ 1 


^\ 


0.0 


y- 


. , \ i 


, , . i / . 


1 
, l\ 


■ 



2b 1.0 




1.5 



?b 1.0 



0.5 



0.0 '. 



+ P 



i i i </ i i 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

v 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

v 







^"=16.971 


H= 


0.5 


2.0 




.. 


i i n i i i • i 


....,..,. 






\ 


\ / 


,-'' 


1.5 




\ 






■ 


et. l.o 








+ p x 


- 


0.5 


T 
1 






- 


0.0 




, 





' 







JT= 


16.971 




H= 


0.0001 




2.0 




1 ' I 




i |i i 


i i i 


|....,V 


" 


1.5 












- 
















- 






" 


Zt. 1.0 






+ P 




- 


0.5 










■ 


0.0 




,-'' 


, i , 


" 


- ; 



0.0 0.5 



1.0 1.5 2.0 2.5 

v 



3.0 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

v 



Figure 1. Zero contours of the functions /i (solid) and /2 (dashed) on the va plane, for different values of the layer depth H. Crossing 
points of the solid and dashed curves represent formal solutions of the dispersion relations. The solution in the limit H — > oo is marked 
by P. 
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Figure 2. Variation of the nondimcnsional growth rate a with layer depth H for M = \2\j2. The solution is analytically fitted by a 
simple shifted tangent hyperbolic (dashed), on which a Gaussian damped sinusoidal is superposed (dotted). 
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Figure 3. Variation of the nondimensional frequency u with layer depth H for J\f = 12\/2. The relation u 2 oc tanh(iJ) (dashed), based 
on an analogy with surface gravity waves, is shown for comparison. 

originating at the interface and the waves reflected from the bottom of the tachocline. To imitate this 
modulation, we superpose an exponentially damped sinusoidal on the tangent hyperbolic to arrive at the 
following fitting function 



a = tanh(H')+0.2sm(l.5H')exp(-0.25H' 2 ), where H' = H - 1.45. 



(23a,b) 



Note, however, that the appearence and details of such interference effects may be sensitive to the details 
of the model setup (i.e. to how sharp the bottom of the tachocline is supposed to be). 

The breakdown of the water wave analogy is even more apparent from figure |3j The variation of the 
frequency with layer depth is clearly quite different for interface dynamo waves than the simple u 2 oc 
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Figure 4. Trajectory of the physical solution sequence in the vo plane for varying H when M = Yl\/2 is kept fixed 

tanh(i^) relation (dashed), valid for surface gravity waves. 

Figure [4] presents the trajectory of the solution in the vfi plane. It is apparent that in the critical case 
v ~ 1.8, at least an order of magnitude higher than the value (2/N) 1 ' 2 suggested by Parker's original 
analysis. 



4 Conclusion 



In this paper we have generalized the textbook interface dynamo model of Parker ( 1993 ) to the case of 
finite tachocline thickness, while otherwise keeping the model as simple as possible. We found that the 
growth rate of the dynamo has a shifted tangent hyperbolic dependence on the layer depth, in agreement 
with the expectation based on an analogy with surface gravity waves. This finding lends some support to 
the claim ( Petrovay]|2007 ) that the combination of a fast tachocline and an interface dynamo can give rise 
to bimodal solutions, possibly explaining the phenomenon of grand minima. Under the assumption that 
the growth rate of the dynamo has a dependence on the tachocline depth similar to what was found in the 
present paper, in that previous work two solutions were found for a nonlinear interface dynamo combined 
with a fast tachocline: a strong field, thin tachocline solution corresponding to normal solar activity and a 
weak field, thick tachocline solution that may potentially represent a grand minimum state of solar acivity. 
Yet the behaviour of the solutions does show some surprising features. The nondimensional frequency 
varies with layer depth in a non-intuitive way and the detailed behaviour of the growth rate — layer depth 
curve displayes a non-monotonic, oscillatory behaviour. This illustrates the richness of phenomena appear- 
ing even in strongly simplified interface dynamo models. A further systematic study of models of this type 
and of more generic models is needed to more fully explore the parameter space and to obtain a deeper 
understanding of the spectrum of possible dynamo solutions in an interface setup. Without such system- 
atic studies no final verdict on the viability of an interface-type solution to the solar dynamo problem is 
possible. 
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